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Abstract 

Optical mushroom-shaped bilhards offer a unique opportunity to isolate and study semiclassical 
modes concentrated on non-dispersive, marginally unstable periodic orbits. Here we show that the 
openness of the cavity to external electromagnetic fields leads to unanticipated consequences for 
the far-field radiation pattern, including directional emission. This is mediated by interactions of 
marginally unstable periodic orbits with chaotic modes. We also show that the semiclassical modes 
are robust against perturbations to the shape of the cavity, despite the lack of structural stability 
of the corresponding classical orbits. 

PACS numbers: 05.45.-a, 05.45.Mt 
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The study of chaotic bilhards has served as a foundation to many areas, ranging from 
nonhnear dynamics and statistical physics to quantum optics. Classically, the phase space of 
a chaotic billiard, defined by its position and momentum variables, is typically divided into 
chaotic and regular regions, which is a property found in generic Hamiltonian systems. How- 
ever, in contrast with smooth dynamical systems, several billiards have been shown to have 
'amilies of marginally unstable periodic orbits (MUPOs) embedded in their chaotic regions 
l|. These orbits, of which bouncing-ball orbits ^ form a particular case, are zero- volume 
structures reminiscent of Kolmogorov- Arnold- Moser islands and dynamically indistinguish- 
able from regular orbits j^. Of the billiards exhibiting this property, no other is attracting 



as much attention as the mushroom-shaped billiards 



Mushroom billiards have the distinctive non-generic feature of exhibiting a single chaotic 



and a single regular region 



J]. These billiards have been used as a model to address the 



impact of MUPOs on chaotic orbits I^, S] and to examine quantum tunneling 0, Q, [s], ^ . 
Yet, the fundamental role of MUPOs in semiclassical dynamics remains largely unexplored. 
The very physical reality of semiclassical MUPO modes remains to be demonstrated given 
that classical MUPOs are not structurally robust against parameter perturbations likely to 



be present in realistic situations 



io|. 



In this Letter, we investigate a family of MUPOs in optical dielectric mushroom cavities. 
"Previous studies have used total internal reflection to identify lasing modes based on stable 
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12l | and unstable periodic orbits [l^ in chaotic cavities, where the latter correspond to 



the so-called scarred modes. In both cases the orbits of interest may not correspond to the 
least leaky modes of the passive system, but give rise to lasing modes via manipulation of 
the gain medium. In contrast, here we show that even in the absence of gain, total internal 
reflection can be used to trap selected MUPOs inside the billiard. All the other orbits leave 
the cavity through refraction, making the MUPO modes the least leaky modes. We explore 
computationally the consequences of this scenario in the semiclassical regime (i.e., short- 
wavelength regime). We show that the reentry of electromagnetic fields into the cavity leads 
to a coupling between MUPOs and chaotic modes, which has no analogue in closed billiards 
and significantly impacts the directionality and intensity of the far-field radiation pattern. 
We also show that the semiclassical MUPO modes are robust against roughness and other 
perturbations to the geometry of the billiard. 

We consider mushroom billiards composed of a circularly shaped hat with radius R and 
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FIG. 1: (Color online) Bottom corner of the hat of the mushroom cavity. The ray represents 
the trajectory of an orbit with reflection angle 9 = 90° — a that hits the undercarriage with 
angle a' < a. When close enough to the corner, the curvature can be ignored and a' = a. The 
inset shows a period-4 MUPO and the entire mushroom cavity with the parameters i?, r and h 
measured from the center of the diameter of the semicircular hat. This point is also the origin of 
the cylindrical and Cartesian coordinates used throughout this Letter, in which z will represent 
the axis perpendicular to the plane of the billiard. 

a foot, square or triangular in shape, that extends to a radial position of r (Fig. [Tj). Regular 
orbits are confined to the hat between the circumferences of radii r and R. MUFOs are also 
confined to the hat but necessarily cross the circle of radius r, while chaotic orbits necessarily 
visit the foot [4] . In the phase space (defined by the coordinates of arc length and reflection 
angle in the circumference of the hat), the regular and chaotic regions are separated by 
the reflection angle 6* = sm~^ (r/R) P, The MUFOs have a reflection angle lying in 
the chaotic region {6 < 6*) but have null Lyapunov exponents, assuring a non- dispersive 
behavior for nearby trajectories. An infinite number of families of MUFOs, characterized 
by different periods and rotation numbers, have been shown to exist for almost all choices 
of r/R . Here we focus on the particular family of period-4 MUFOs (Fig. [1]), which 
exist for r/R> l/\/2 and have a reflection angle 9p = 45°. 

These MUFOs can be trapped and isolated in an open mushroom cavity with certain 
refractive index n by exploiting the critical incidence angle 6c = sin~^(l/n) below which 
light escapes the cavity. Isolating the target MUFOs requires that all chaotic and regular 
modes leave the cavity. The former is easily accomplished since chaotic orbits will eventually 
have a reflection angle smaller than the critical angle, as in previous applications to select 
stable modes in chaotic resonators \v^. The latter, on the other hand, can be accomplished 
by forcing the regular orbits to escape through the "undercarriage" (bottom of the hat), as 
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demonstrated next. 

Figure [1] shows the trajectory of an orbit. If the reflection is close enough to the corner, 
the curvature of the hat can be ignored and a' = a. The MUPOs and the refractive index 
are chosen such that the reflection angle satisfies > 6c and a'^ = 90° — 6p> 9c. We also 
want a cavity in which a^gguiar < ^c- The smallest reflection angle for any regular orbit in 
the mushroom cavity is 9* . Therefore, our conditions become 90° — 9* < 9c < 90° — 9p and 
9c < 9p. 

We shall consider two commonly used materials in microcavity fabrication: polymer and 
semiconductor (Table [11). We first focus on polymer, which has refractive index n = 1.5 
and critical angle 9c ~ 41.8°. We assure 90° — 9* < 9c by taking r/R = 0.75, which yields 
90° — 61* ^ 41.4°. Thus, since 9p = 45° > 9c, our conditions are met and the period-4 MUPOs 
will be isolated in the cavity. 

TABLE I: Parameters of the materials, period-4 MUPOs, cavities, and FDTD grid spacing con- 
sidered. 





n 9c Op 


9* 


r/R 


h/R 


A/(nAx) 


polymer 


1.5 41.8° 45° 


48.6° 


0.75 


0.70 


20 


semiconductor 


3.3 17.6° 45° 


46.9° 


0.73 


0.68 


20 



The above argument assumes we are in the classical regime and is expected to be approx- 
imately valid in the semiclassical regime considered here, which is characterized by nkR ^ 1 
and hence small wavelengths A = 27r/k compared to the size R of the cavity. We set the 
parameters of the cavity to be i? = 10/im, r = 7.5/im, and h = 7^m. We consider TE 
polarization [E^ = 0) and employ a parallelized version of the finite-difference time-domain 
(FDTD) method iJ] in 2D with grid spacing of Ax = Ay = 1.36 nm for a square grid in 
the plane of the cavity. The cavity is excited uniformly with a sinusoidal-Gaussian pulse 
and the simulation run until only the eigenmodes with the highest Q = uj/'j remain, where 
u is the angular frequency and 7 is the energy decay rate of each eigenmode. 

The eigenmodes are determined by a Fourier transformation of the energy in the cavity. 
A wide range of excitation wavelengths were simulated before arriving at the highest Q 
profile for A = 41.0 nm. The high Q modes are limited to such small wavelengths most 
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FIG. 2: (Color online) Spatial intensity for the period-4 MUPO mode with A = 41.0 nm and 

Q = 13 000 in a polymer microcavity. The intensity was averaged over one oscillation period. 



likely because of the distance the MUPO modes are restricted to traveling in the corner of 
the mushroom hat. 

Because we work with the TE polarization, our far- field analysis deals with Hz- The 
field outside the cavity as a function of the angle at radial distance p from the center of 
the diameter of the hat behaves as Hz{4>) = arp,Hmikp)e'^^'^, where Hm is the Hankel 
function of the first kind with azimuthal order m 15|. The coefficients am are found by 
applying 27iamHm{kp) = A0 ^ ■ i7^(0j)e~™'^J , where (pj is the angle discretized according 



to a resolution A(f). The field is then evaluated at a far-field distance p/, determined by 
the Fraunhofer condition pf ^ -R^/A, to obtain the far-field pattern In general 

Hz{4>j) is complex, but FDTD only calculates the real component of the field. We obtain 
the imaginary component from Faraday's law ImHz = {—ujpo)^^^ x E, where E is real. 

Figure [2] color-codes the field intensity for a polymer cavity with a triangular foot, indicat- 
ing that the highest Q mode is in this case a single period-4 MUPO mode with A = 41.0 nm. 
A fit to the mono-exponential curve of energy versus time yields a decay rate 7 that re- 
sults in Q = 13 000. This eigenmode was isolated by narrowing the excitation width to 
6X = 0.005 nm. In closed billiards, the shape of the foot does not play a crucial role. In the 
open cavities considered here, however, the shape of the foot can significantly infiuence the 
dynamics of the modes. This is particularly so for polymer cavities since their modes are 
prone to leak. For example, transitioning from a triangular foot to a square foot results in 
the MUPOs no longer being the highest Q modes IGj: new high Q modes are established via 
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FIG. 3: (Color online) (a) Far-field radiation pattern |ff/(</>j)p for the polymer microcavity and 
high Q period-4 MUPO mode shown in Fig. El (b) Magnification of |i?/((/>j)p for 0° < < 30°. 
The far-field pattern is strongly bidirectional. The angular resolution is A(j) = 0.008°. 

leakage from the undercarriage and subsequent penetration into the foot. This is expected 
to be relevant to optical experiments because the highest Q mode typically determines the 
first lasing mode. More strikingly, when the MUPO remains the highest Q mode, as in the 
triangular-foot example of Fig. [2l field reentry can lead to directional emission. 

Figure [3] shows the far-field pattern generated by the period-4 MUPO mode shown in 
Fig. [21 The emission pattern is strongly bidirectional, with peaks at approximately 10° and 
170°. The percentage U of total emission within [8.5°, 11°]U[169°, 171.5°], is U ^ 25% which 
is 18 times the corresponding value for isotropic emission. From our prior discussion of 
internal reflection, much of the field leakage should be in the form of an evanescent field. 
When the field hits the top of the hat, the evanescent field outside the cavity is diffracted at 
various angles due to the curvature. A similar event takes place as the evanescent field along 
the undercarriage diffracts at the corner of the hat. In both cases, there is no significant 
contribution to the primary peaks of the far-field pattern. 

The directional emission is actually caused by refraction through the undercarriage and 
subsequent penetration into the foot. As suggested in Fig. [2] and confirmed in Fig. HI this 
is so because part of the MUPO's field strikes the boundary further up the curvature of the 
hat and hits the undercarriage with an incidence angle a' < Oc- The field then refracts out 
of the cavity towards the foot. Once the field enters the foot, it is refracted upwards to the 
boundary of the hat and immediately escapes. Figure ll](a) shows the corresponding classical 
trajectory with d marking the ray's incident location on the undercarriage, (3 the upward 
angle of refraction once the ray reenters the cavity, and (p the emission angle. Figure IH^b) 
shows the values of a' and d that traced along this trajectory satisfy 8.5° < < 11°. Only 
incidence angles close to the critical angle — with a' = 38° being the minimum value found — 
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FIG. 4: (Color online) (a) Calculated ray diagram indicating the point of incidence on the under- 
carriage d, reentry angle /3, and emission angle in a mushroom cavity, (b) Values of a' and d 
that result in emission angles 8.5° < (j) < 11° (shaded area). The corresponding values of /3 (black 
curve) agree with the results shown in Fig. [2j 

result in emission angles agreeing with the far-field pattern shown in Fig. [3l Moreover, the 
smallest incidence angles only result in agreement when d > 8.8/im, which is consistent with 
Fig. [2] in that the mode is concentrated in the corner of the hat. A measurement of the 
reentry angle from the Fig. [2] data gives (3 ~ 13.5° (with some spreading), which is also 
consistent with the calculation of (3 in Fig. Hl^b). 

We now turn to semiconductor microcavities, which are computationally less demanding 
and less leaky and, as shown below, have high Q modes in the optical range and can be 
used to further validate the mechanisms found in the polymer cavity. With a refractive 
index n = 3.3, the semiconductor has a smaller critical angle of 6c ~ 17.6°. We consider 
this case for R = lO/im, h = 6.8/im, and r = 7.3/im, so the period-4 MUPOs are located 
just inside the chaotic region; we use a larger grid spacing of Aa; = 8.3 nm and continue to 
focus on TE polarization (Table [T]). One of the conditions for isolating the MUPO modes 
is no longer satisfied, namely 90° — 6* ^ 9c, rendering some of the regular orbits to never 
reach the undercarriage with incidence angle smaller than 6c. Nevertheless, because their 
incidence angles are farther apart from the critical angle by at least 6* — 6p = 1.9°, the 
period-4 MUPOs remain the highest Q modes. 

Indeed, with an excitation wavelength of A = 580 ± 25 nm, the highest Q eigenmodes are 
excited in the range 500 nm < A < 600 nm and are all found to be period-4 MUPO modes 
irrespective of the foot shape. A harmonic inversion software [l7| was used to estimate the 
Q values. The highest Q modes were selected and simulated individually for a more accurate 
estimate of Q. A larger wavelength MUPO mode at A = 598 nm has Q = 34 000 while a 
smaller wavelength MUPO mode at A = 526 nm has Q = 160 000. Because resolution and 
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FIG. 5: (Color online) (a) Fourier transformation for the period-4 MUPO modes excited by A = 
550.92 lb 0.5 nm in a semiconductor microcavity. There are two eigenmodes within this excitation 
range, at wavelengths A = 550.61 nm and A = 550.91 nm. The Q values for these modes are 39000 
and 140 000, respectively. (b,c) Far-field pattern \Hl{(f)j)\^ corresponding to (b) A = 550.61 nm 
and (c) A = 550.91 nm. 

simulation runtime are an issue for smaller wavelengths, the two eigenmodes identified in 
Fig. (S^a), at A = 550.91 nm with Q = 140 000 and at A = 550.61 nm with Q = 39 000, are 
chosen for further investigation. These two modes are representatives of the typical behavior 
found for the semiconductor cavity. As shown in Figs. [5]^b)-(c), the corresponding far- field 
pattern is close to isotropic for A = 550.61 nm and slightly bidirectional for A = 550.91 nm. 

Figure [6] shows the spatial intensity for these modes. The low Q mode at A = 550.61 nm 
is concentrated close to the corner of the hat [Fig. [6](a)]. The evanescent field strongly 
diffracts at the corner resulting in a nearly isotropic far-field pattern. This case is thus 
similar to the one in Fig. [2] and provides further evidence that the directionality in the 
polymer microcavity originates from field reentry and subsequent refraction as opposed to 
direct refraction or diffraction of evanescent waves. The high Q mode at A = 550.91 nm, 
on the other hand, is concentrated much closer to the intersection of the undercarriage and 
the foot, and hence has less leakage from diffraction at the corner [Fig. [6]^b)]. The slight 
bidirectionality of the high Q mode observed in Fig. [5]^c) is caused by scattering of the field 
near the intersection. 

The robustness of the high Q period-4 MUPO mode shown in Fig. EJ^b) was tested 
by rounding the corners of the cavity with a radius of curvature Re and by introducing 
roughness to the radial boundary of the hat and undercarriage. The latter is defined by 
a standard deviation A and modulation period D. As the perturbation increases, either 
through rounded corners or roughness in the boundary, Fourier transformations reveal that 
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FIG. 6: (Color online) (a,b) Spatial intensity \Hz\'^ for the period-4 MUPO modes at (a) A = 
550.61 nm, Q = 39 000 and (b) A = 550.91 nm, Q = 140 000 in a semiconductor microcavity. The 
intensity was averaged over one oscillation period. Different times are selected to yield the same 
maximum intensity for both plots. 



the wavelength of the MUPO mode changes little. The largest roughness for which the 
MUPO mode is not destroyed is for A = 9nm and D = 40A. Even in this case, the 
resulting wavelength shift was just 0.01%. Similar or even more pronounced robustness was 
observed for the directionality of the far-field pattern of the highest Q mode of the polymer 
microcavity [16,]. However, the Q values proved quite sensitive to the perturbations. For 
the case Fig. Mjo), the Q value reduces from 140 000 to 73 000 for Re = 300 nm and to 
Q = 3 200 for A = 9 nm and D = 40A. Nevertheless, the spatial intensity pattern of the 
MUPO mode is fairly well preserved in all cases considered. Indeed, unlike the classical case, 
where a MUPO either exists or does not exist, the semiclassical MUPO mode changes with 
the perturbations but mimics the unperturbed MUPOs for the entire lifetime of the mode 
in the cavity. 

MUPOs are inherently less dispersive than unstable periodic orbits and, as shown here, 
underlie modes that can be isolated and trapped in open optical cavities. The openness and 
non-convexity of the cavity, and consequent reentry of refracted electromagnetic fields, can 
lead to directional emission, fundamentally altering the far-field radiation pattern. The re- 
cently demonstrated prevalence of MUPOs in many billiards [l| indicates that these findings 
are most likely not limited to the mushroom billiards considered in this Letter. 
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